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Power  Spectra  of  Internal  Gravity  Waves 


1.  INTRODUCTION 

This  report  answers  a  question  implied  In  three  earlier  publications  (Dewan,  1979;  Dewan 
and  Good,  1984,  1986).  The  question  is  "Can  the  so-called  universal  spectral  model  of  Garrett 
and  Munk  (1972,  1975,  1979)  be  explained  by  physical  similitude  arguments?"  The  apparent 
success  of  the  similitude  argument  In  Dewan  and  Good  (1984,  1986)  suggests  that  it.  Indeed,  can 
be  done.  The  original  models  of  Garrett  and  Munk  (GM)  cited  above  were  based  on  a 
combination  of  empirical  observations  of  oceanic  internal  gravity  waves  plus  the  general 
properties  of  these  waves  that  can  be  deduced  from  the  equations  of  motion  such  as  dispersion 
relations  and  so-called  "wave  functions".  Subsequently,  VanZandt  (1982)  showed  that  the 
universal  spectral  model  also  applies  to  the  internal  waves  in  the  atmosphere.  The  term 
"universal"  in  both  contexts,  however,  must  be  understood  as  somewhat  of  an  exaggeration. 

For  example,  the  amplitude  of  the  spectra  fall  "mostly  within  a  factor  of  2"*  (GM,  1972,  GM, 
1975)  and  the  log-log  slopes  of  certain  wave  number  spectra  are  between  5/3  and  2  (GM.  1972) 
or  2.5  (GM,  1975)  or,  in  the  case  of  the  atmosphere,  nearly  3  (Dewan  et  al,  1984). 

In  an  earlier  paper  (Dewan,  1979)  it  was  suggested  that  the  GM  spectral  model  might  be 
applicable  to  the  atmosphere;  but,  the  situation  studied  in  that  paper  involved  exceptionally 

(Received  for  publication,  17  September  1990) 

In  Section  6  it  will  be  mentioned  that,  in  the  atmosphere,  a  change  of  altitude  can  cause 
the  amplitude  of  the  temporal  PSD  to  change  by  a  factor  of  200. 
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vigorous  waves  In  the  stratosphere.  No  attention  was  given  there  to  the  "universal"  amplitude 
property:  but  it  was  suggested  that  the  waves  Involved  were  Interacting  sufficiently  to  provide  a 
"cascade  effect"  to  explain  spectral  slope.  An  analogy  to  turbulent  cascades  was  made  there. 

Historically  speaking,  the  turbulent  cascade  was  first  described  In  a  poem  by  Richardson 
(1922):  "Big  whirls  have  little  whirls  that  feed  on  their  velocity,  and  little  whirls  have  lesser 
whirls,  and  so  on  to  viscosity  (in  the  molecular  sense).”  In  Dewan  (1979)  a  detailed  argument 
was  presented  for  the  case  that  a  similar  thing  could  be  occurring  for  atmospheric  gravity 
waves  and  that  "Big  waves  have  little  waves  that  feed  on  deformation,  and  little  waves  have 
lesser  waves  to  turbulent  dissipation  (In  the  eddy  sense)."  Kolmogorov  deduced  his  Inertial 
range  spectrum  for  velocity  fluctuations,  <j>  =  Cj  e2^3  k  -5/3  (e  is  the  dissipation  rate,  k  the  wave 
number  and  cta  universal  constant  of  order  unity)  by  means  of  a  similitude  argument  based  on 
the  Richardson  poem.  In  Dewan  (1979)  4>  =  c2  e2^3  1^  ~5^3  (kx  Is  horizontal  wave  number  and 
c2  *Cj)  was  derived  In  exactly  the  same  manner  via  the  second  "poem".  In  this  report  1  wish  to 
extend  the  argument  In  Dewan  (1979)  in  order  to  attempt  to  find  a  physical  explanation  of  the 
remainder  of  the  GM  spectral  model. 

In  apparent  contradiction  to  the  above  approach,  Dewan  and  Good  (1984,  1986)  explained 
the  universal  spectrum,  Tu  (k^)  =  c3  N2  k^3  (where  S'  (kz)  is  the  PSD  for  horizontal  velocity 
fluctuations,  N  Is  the  buoyancy  frequency  and  kz  the  vertical  wave  number)  which  Is  observed 
in  the  stratosphere  (Dewan  et  al,  1984),  and  troposphere  (Daniels,  1982,  and  Endlich  and 
Singleton,  1969)  by  means  of  an  argument  based  on  "saturation"  rather  than  on  "cascade” 
effects.  Again,  a  similitude  argument  was  basic  to  this  approach  in  spite  of  the  fact  that 
several  other,  more  detailed  arguments  were  also  presented.  This  work  was  subsequently 
augmented  by  Smith  et  al  (1987)  who,  among  other  things,  deduced  the  value  c3  above  and, 
most  importantly,  predicted  the  altitude  dependence  of  the  dominant  vertical  wave  number 
which  they  called  m.. 

In  subsequent  years  there  have  been  further  tests  of  the  "saturation  hypothesis"  (in  the 

form  of  the  S^  (kz)  =  N2  k^3  dependence).  For  example  see  Fritts  et  al  (1988),  Gardner  et  al 
(1989),  Gardner  and  Voelz  (1987),  Gardner  (1989),  Hass  and  Meyer  (1987)  and  Wu  and  Wlddel 
(1989).  More  recently  Tsuda  et  al  (1989)  wrote:  'The  fact  that  the  argument  of  the  saturated 
gravity  wave  spectral  model  of  Dewan  and  Good  (1986)  and  Smith  et  al  (1987)  improves  as  the 
observed  spectra  are  Improved  suggests  that  the  model  is  essentially  correct,  in  spite  of  the 
heuristic  nature  of  some  of  Its  assumptions." 

Most  recently,  Hines  (1990)^  re-examined  these  "heuristic  assumptions"  and  seems  to  have 
put  the  entire  program  on  a  more  extended  physical  footing.  In  addition,  he  also  raised  the 
question  that  the  saturation  picture  may  give  the  correct  answer  for  the  wrong  reason.  While 
this  seems  not  to  have  disturbed  the  similitude  argument  in  Dewan  and  Good  (1984,  1986), 
Hines'  new  physical  picture  is  very  stimulating  to  further  progress  with  this  problem.  The 


1  A  detailed  response  to  the  criticisms  of  Hines  (1990)  is  postponed  to  a  future  publication. 
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reader  is  encouraged  to  examine  Hine's  paper  for  the  details;  but,  his  approach  involves 
strongly  interacting  waves  and  the  "compliant  wave  number"  regime  first  discussed  by  Munk 
(1981).  In  this  regime  the  small  waves  (m  >  mj  are  doppler  shifted  and  strongly  modified  by 

critical  layer  processes  due  to  the  convective  effects  of  the  large  waves.  In  a  word,  Hines’  new 
picture  has  made  one  thing  very  clear  for  the  present  purpose,  and  this  is  that  the  success  of 

the  HP  (kz)  =  c3  N2  k^3  spectrum  no  longer  argues  against  cascade  processes  which  is  ours.t  In 
fact  the  opposite  is  true;  that  is  to  say,  the  cascade  effects  in  waves,  as  suggested  in  Dewan 
(1979),  appear  more  valid  when  viewed  from  the  perspective  of  Hines'  "wave-interaction" 
suggestion.  Usually  one  regards  "saturation"  as  implying  that,  as  waves  grow,  they  eventually 
reach  a  "saturated  amplitude"  above  which  they  "dump  energy"  directly  from  their  large  scales 
into  small  turbulent  scales.  In  contrast,  Hines'  picture  of  strongly  interacting  waves  is 
consistent  with  cascade  processes  of  energy  flowing  from  large  to  small  scale  via  intermediate 
scales.  A  second  benefit  here  from  the  Hines  picture  is  its  agreement  with  the  long  discussed 
and  disturbing  fact  that  gravity  waves  are  not  the  weakly  interacting  waves  that  theorists 
would  prefer  them  to  be  for  calculational  convenience.  Rather  they  are  unavoidably  strongly 
interacting.  The  latter  point  has  been  emphasized  in  the  oceanic  internal  wave  context  by 
Holloway  (1980,  1981). 

For  completeness  it  should  be  mentioned  that  certain  previous  authors  have  raised  the 
question  of  whether  or  not  the  GM  model  or  the  saturated  spectrum  is  related  to  buoyancy 
range  (BR)  turbulence.  Weinstock  (1985),  for  example,  has  suggested  that  the  saturated 
spectrum  may  be  the  (BR)  turbulence  as  described  by  Lumley  ft  (1964),  which  also  has  N2k~3 
dependence.  This  view  was  also  suggested  by  Holloway  (1983)  and  Gargett  (1981).  Also,  Sidi  et 
al  (1987)  showed,  on  the  basis  of  saturation,  that  the  temporal  part  of  the  GM  spectrum  would 
have  a  of1  dependence  and  he  considered  this  as  a  form  of  BR.  In  what  follows  I  will,  in 
contrast,  get  a  co-2  dependence,  but  in  general,  I  am  at  this  point  somewhat  in  overall 
agreement  with  the  general  idea  that  the  GM  spectrum  may  be  of  "BR  turbulence".  On  the  other 
hand,  the  word  "turbulence"  in  this  context  is  very  misleading  in  the  sense  that  one  usually 
assumes,  by  definition,  that  turbulence  does  not  obey  dispersion  relations.  As  discussed  in 
Dewan  (1985),  there  may  exist  fluctuations  that  are  between  the  extremes  of  turbulence  and 
waves  which  continue  to  retain  dispersion  relations  and  limited  propagation  properties,  and 
yet  have  very  active  nonlinear  mode  interactions.  Perhaps  this  is  what  physically  happens  in 
BR  "turbulence"  as  well  as  in  the  "universal  spectrum  range  of  gravity  waves";  and.  it  would 
explain  why  similitude  arguments  based  on  cascades  are  successful  (if  indeed  the  main 
predictions  I  make  here  are  found  to  be  valid  in  future  experiments). 

In  order  to  relate  the  GM  formalism  to  similitude  arguments  I  will  first  give  a  brief  review 
of  their  formalism  in  Section  2  and  then  present  the  derivation  of  the  temporal  part  of  their 
formalism  in  Section  3.  In  Section  4  a  number  of  testable  PSD's  are  derived  which  are  then 


t  On  the  other  hand,  since  it  refers  to  interaction  between  a  single  large  scale  region  and  a 
single  small  scale  region,  it  does  not  necessarily  imply  that  a  cascade  is  involved. 

tr  Compare  Phillips,  p.  293  (1980)  whose  discussion  is  highly  to  be  recommended  especially 
in  conjunction  with  what  he  says  on  p.  232  regarding  "promiscuous"  mode  interaction. 
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collected  in  Section  5  where  predictions  for  future  tests  are  described.  In  Section  6  a  limited 
comparison  between  theory  and  experiment  is  made  which  demonstrates  that  the  predictions 
are  not  unreasonable  in  the  sense  that  they  appear  to  be  of  the  right  order  of  magnitude. 

Finally,  still  further  experimental  support  is  obtained  by  comparing  the  altitude  dependence  of 
a  temporal  PSD  of  velocity  fluctuations  with  the  altitude  dependence  of  dissipation  rate,  e, 
since  the  theory  says  that  there  is  a  direct  relationship  between  them. 


2.  A  BRIEF  REVIEW  OF  THE  GARRETT-MUNK  (GM)  FORMALISM 


The  GM  (1972,  1975)*  formalism  for  the  spectral  model  of  internal  oceanic  gravity  waves 
begins  with  the  spectrum  for  scalar  energy  per  unit  of  mass  which  they  designate  by  E  (kz,  to) 
where  kz  is  the  vertical  wave  number  and  to  is  the  radial  frequency  (kx,  the  horizontal  wave 
number  is  also  used).  They  represent  this  spectrum  as  being  factorable: 


E  (kz,  to)  s  E'  B(to) 


(2.1) 


where  k  represents  the  vertical  wave  number  for  the  scale  where  the  energy  begins  to  decrease 
z« 

rapidly  with  kz.  It  also  is  the  "dominant"  or  "characteristic"  wave  number  in  the  sense  that,  as 
Smith  et  al  (1987)  have  shown,  "most  of  the  energy  lies  within  a  factor  of  3  of  the  mode  at 

(k  /21/3)."  The  Ms  a  non-dimensional  wave  number: 

z- 


*-(VkJor(Vk«.} 


(2.2) 


Also  A(X)*(t-  1)  (1  +X.)  \ 

In  GM  (1975)  the  value  t  =  2.5  was  chosen.  In  their  formalism  non-dimensional  quantities 
were  used  by  making  use  of  normalization  with  representative  scales  of  time  and  distance. 

For  example,  all  frequencies  were  normalized  by  N  where  N  is  the  "extrapolated  Vaisal  a 

frequency"  in  their  depth  dependent  buoyancy  frequency  profile  given  by  n  =  Ne"^  h  where  n  is 
the  dimensional  local  frequency,  y  the  depth  and  b  the  characteristic  length  for  this  depth 
dependence. 


*  More  recent  versions  appeared  in  GM  (1979)  and  Munk  (1981). 
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The  "hats"  over  the  letters  were  used  by  GM  to  designate  dimensionality^  Wave  numbers 

/N  -J 

were  normalized  by  M  ■  b  2  n  In  the  following  I  will  try  to  avoid  use  of  their  non- 
dimensional  approach  because  (a)  in  the  case  of  the  atmosphere  it  Is  known  that  this  would 

✓v 

give  an  Incorrect  scaling  with  respect  to  "N".  and  (b)  the  dimensional  approach  is  more 
convenient  for  the  experimental  comparisons  that  will  be  made. 

Both  B(w)  and  A(X)  are  defined  so  they  integrate  to  unity.  This  will  be  equally  valid  even 
when  dimensional  terms  are  used  because  the  Integrals  themselves  are  dimensionless.  Thus 
we  have 


B(o>)  dro  =  1;  and 


f 


A(X)  dX  =  1 


(2.3) 


where  cot  Is  the  inertial  frequency.  In  order  to  obtain  the  spectra  ("power  spectral  densities”  to 
be  more  exact,  henceforth  denoted  by  "PSD's")  of  velocity  components,  GM  use  what  they  call 
"wave  functions"  in  a  manner  to  be  shown  below.  In  addition,  they  Integrate  over  X  or  w  in 
order  to  obtain  the  one  dimensional  PSD's.  For  example,  in  order  to  obtain  the  PSD  with 
respect  to  ltz  of  the  horizontal  velocity  fluctuations,  u,  they  use 


¥ 

U 


(\)  = 


_|N  E'B((o)A(k/kzJ 


~2 

U  (o))do> 


(2.4) 


where 


~2 

U  (to) 


n 


(«•♦■?) 

2 

<0 


(2.5) 


~2 

and  where  U  (o>)  is  the  wave  function  for  horizontal  velocity,  again,  co(  is  the  inertial 
frequency,  E'  is  a  dimensionless  energy  which  has  been  empirically  determined,  and  n  is  a 
dimensionless  parameter  relating  to  buoyancy  frequency  which  I  shall  set  =  1  for  the 
atmospheric  model.  The  correct  dependence  on  N  Is,  from  Dewan  and  Good  (1984,  1986), 


*  In  this  paper  I  shall  use  this  notation  exclusively  to  refer  to  the  GM  model.  The  new  model 
being  presented  assumes  that  all  quantities  have  dimensions. 
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a  N2kz  for  the  universal  part  (l.e.  high  wave  number)  of  the  spectrum.  Smith  et  al  (1987) 
brought  this  Into  the  context  of  the  GM  formalism  by  writing  It  as 


¥  (k  )  = 

u  V  z)  3 


z. 


f  k  \ 

1  +  — £ 
k 

v  z»y 


-3 


(2.6)t 


Notice  that  It  Is  here  that  the  dependence  on  N  differs  from  the  GM  model,  but  Is  supported  by 
atmospheric  observations  (Dewan  and  Good,  1986).  Smith  et  al  (1987)  found  a  value  for  a,  as 
has  been  mentioned,  and  It  is  a  =  1/6  (compare  Hines,  1990). 

The  vertical  velocity,  v,  wavefunction  is  given  by  GM  (1972,  1975)  as  (n  =  1) 


_ 2 

V  (co)  = 


N2 


(2.7) 


Thus 


'  B<«>>  A(Vk0 


2 

(co)do). 


Z. 


(2.8) 


Similar  relations  hold  for  the  kx  type  PSD’s  where  in  that  case  one  uses  X  =  k^/k^ 
The  (o  dependence,  for  example,  Is  given  by 


'*(«>)  =  U  (co)E'B(co) 


f 


A  (X.)  d^. 


(2.9) 


=  U  (w)E'  B(co) 


from  Eq.  (2.3). 


t  In  actual  fact.  Smith  et  al  (1987)  used  a  slight  modification  of  this  for  convenience. 
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In  order  to  obtain  PSD's  for  displacements  one  employs  the  relevant  velocity  component 
wave  function  and  divides  It  by  co2.  For  example,  the  PSD  for  vertical  displacement,  as  a 
function  of  co,  is  given  by 


~2 


0) 


(2.10) 


and  so  on. 

It  should  be  mentioned  that  Garrett  and  Munk  (GM)  (1972,  1975)  obtained  these  wave 
functions  In  a  complicated  way;  but,  they  can  be  obtained  directly  from  the  polarization 
relations  when  suitable  approximations  are  made. 

The  values  of  E  and  B(co)  In  the  GM  models  are 


/s 


(2.11)* 


and 


E  =  2rt x  10  5N 


/M 


(2.12) 


(l.e.  E'  =  2n  x  10"5). 

In  the  next  section  I  will  deduce  E  B(co)  from  a  similitude  argument.  The  above  model  is 
strictly  an  "empirical  fit".  In  both  cases,  It  should  be  pointed  out,  a  "non-equatorial"  model  is 
assumed.  This  is  because  something  unphysical  happens  when  one  sets  cot  =  0. 


t  Here  co,  and  co  are  dimensional. 
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3.  A  SIMILITUDE  ARGUMENT  FOR  E  B(co)tt 


As  was  mentioned  previously,  a  rationale  has  been  given  In  Dewan  (1979)  for  an 

—5/3 

atmospheric  gravity  wave  cascade  and  that  this  leads  to  the  observed  kx  dependence.  In  the 
present  section  the  same  strong  mode-interaction  cascade  picture  will  be  used  to  derive  E  B(<o). 
The  cascade  assumption  made  here  will  also  entail  an  additional  premise,  that  is,  the  wave 
cascade  Is  In  steady  state.  In  other  words,  I  shall  assume,  as  in  Dewan  (1979)  and  as  in  the 
case  of  inertial  range  turbulence,  that  the  input  and  output  of  the  cascade  are  both  equal  to  the 
dissipation  rate,  e,  and  remain  approximately  constant  over  time. 

These  assumptions  make  It  easy  to  obtain  E  B(co),  for,  the  latter  energy  spectrum  can  only 
depend  upon  two  parameters,  namely  co  and  e.  These  have  the  dimensions  of  T"1  and 
L2  T~3  respectively.  The  dimensions  of  E  B(co)  are  of  velocity  variance  per  frequency  band¬ 
width  or  L2!^1.  Thus 


L  ,  a  b 
—  =  a  e  co 
T 


(3.1) 


For  equality,  a  *  1 ,  b  =  -2  thus 


E  B(co)  =  a' 

co 


QED. 


(3.2) 


where  a'  Is  a  universal  constant  of  order  unity.  (Here  B(co)  =  (co,  /  co2)  and 
E  *  (a'e  /  co,).] 

It  is  not  surprising  that  this  is  Identical  to  the  result  found  in  Tennekes  and  Lumley  (1972) 
for  the  inertial  range  frequency  spectrum.  More  interesting  is  the  fact  that,  In  the  turbulence 
case,  Eq.  (3.2)  does  not  agree  with  experiment.  In  view  of  this,  Tennekes  (1975)  subsequently 
did  a  more  careful  analysis  of  isotropic  Inertial  range  turbulence  In  which  he  took  into 
account  the  advective  doppler  effects  of  large  eddies  upon  the  measured  frequencies  of  the 
small  eddies  when  the  latter  are  measured  in  the  Eulerian  frame.  He  obtained 


<K©)  -  P  e 


2/3q2/3o)-5/3 


(3.3) 


ft  The  hat  notation  is  dropped. 
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(where  1  /2  q2  is  the  mean  kinetic  energy  per  unit  mass}  which  does  agree  with  experiment.  As 
will  be  discussed  below,  in  the  case  of  oceanic  gravity  waves  the  ©-2  dependence  seems  to  be 
correct,  whereas,  from  VanZandt  (1982),  in  the  case  of  the  atmosphere  there  is  evidence  for 
(o~5/3.  xhe  Tennekes  argument  offers  a  way  to  explain  the  difference,  and  it  would  entail 
verifiable  predictions  regarding  q.  For  the  details  see  the  Appendix. 

In  Dewan  and  Good  (1986)  physical  interpretations  of  the  similitude  results  were  given.  As 
was  mentioned,  these  interpretations  have  been  questioned  by  Hines  (1990).  Similar  physical 
pictures  for  the  cascade  are  available.  For  example,  one  can  imagine  that  each  band-width 
around  a  given  frequency,  to,  has  an  associated  velocity  variance  which  I  will  call  v2!©).  The 
latter  is  related  to  e  and  t(co),  the  "mode  interaction  time"  of  the  ©  mode.  In  other  words,  t(©) 
is  the  time  needed  for  the  "©-band"  in  question  to  give  up  a  significant  fraction  of  its  "kinetic 
energy  per  unit  mass",  v2!©).  We  are  here  assuming  strong  interactions.  In  this  case  t  begins  to 
become  somewhat  "comparable"  to  the  period  of  the  mode  in  question  (compare  Dewan,  1985). 

I  therefore  let  x  ~  ©_1.  Thus  v2!©)  would  scale  as  e/©.  Next  I  make  an  assumption  similar  to 
one  found  in  Dewan  and  Good  (1986)  and  let  the  band-width,  A©,  scale  as  ©.  Then,  1  find  from 
this  that 


(const.) 

© 


(3.4) 


as  in  Eq.  (3.2).  But  this  argument  of  A©  ~  ©  is  exactly  the  sort  of  thing  Hines  (1990)  finds 
inadequate  because,  left  unanswered  are  some  questions:  (a)  why  is  A©  ~  ©?  If  the  answer  is 
that  you  have  wave-packets  of  frequency,  say,  ©,  and  that  they  all  contain  a  comparable 
number  of  oscillations,  then  the  next  questions  are  (b)  why  should  that  be  true,  and  (c)  why 
should  all  the  wave-packet  scales  get  "equal  representation"  in  time?  To  avoid  such  questions 
we  must,  for  now  at  least,  let  the  similitude  argument  stand  on  its  own  feet.  This  is  regrettable 

because  it  will  be  shown  below  that  in  this  model,  E  =  — -  (which  makes  a  lot  of  physical  sense 

© 

i 

in  my  cascade  picture  since  most  of  the  energy  is  at  frequency  ©,  and  hence  v2^)  ~  e/©, ) 
represents  most  of  the  ”E"  to  be  found.  Is  it  really  impossible  to  answer  the  embarrassing 
questions  raised  by  Hines? 

A  way  out  of  this  dilemma  was  suggested,  in  private  communication,  by  D.  Fritts.  His 
suggestion  is  that  a  "dynamical  similarity"  exists  at  all  relevant  scales  "like  fractals".  Fritts' 
suggestion,  it  should  be  pointed  out,  is  totally  consistent  with  the  cascade  picture  espoused  in 
the  present  paper.  In  fact,  such  fractal  scaling  should  be  directly  caused  by  the  cascade 
process.  Put  another  way,  the  relationship  between  self-similar  cascades  in  turbulence  and 
self-similar  fractals  is  well  known.  See  for  example  Mandelbrot  (1975  and  1983).  The  point 
here  is  that  it  is  postulated  that  waves  too  can  have  self- similar  cascades  (Dewan,  1979  and  the 
discussion  above)  and  this  gives  rise  to  the  self-similarity  conditions  A©  ~  ©  and  Ak  ~  k,  and 
the  "equal  representation"  assumptions  which  were  previously  left  unexplained  in  Dewan  and 
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Good  (1986)  and  Smith  et  al  (1987).  The  cascade  hypothesis,  therefore,  answers  at  least  the 
main  questions  and  objections  raised  by  Hines  (1990). 

Experiments  suggested  below  involve  simultaneous  measurements  of  PSD’s  and  e  as  well  as 
k^,  etc.  If  the  predictions  of  the  cascade  model  hold  up,  then  the  present  "theory"  will  have 
been  "verified".  It  would  be  desirable,  then,  to  go  beyond  similitude  arguments;  however,  as  in 
the  case  of  inertial  range  turbulence  it  is  not  entirely  clear  how  this  could  be  done  rigorously. 
Dimensional  analysis,  or  similitude,  can  sometimes  be  surprisingly  powerful  and,  at  the  same 
time,  very  simple.  It  would  constitute  a  grave  error  to  mistake  this  simplicity  for  lack  of  rigor 
or  substance  in  my  view.  The  Inertial  range  turbulence  model  of  Kolmogorov,  which  is  based 
eiitirely  on  similitude,  is  the  most  rigorous  and  experimentally  valid  model  for  this  type  of 
turbulence  yet  published.  It  has  stood  the  test  of  time  for  one-half  century. 


4.  THE  STRONG-MODE-INTERACTION  MODEL  FOR  THE  "UNIVERSAL  SPECTRUM" 


A  comprehensive  PSD  model  can  be  obtained  from  the  results  of  Section  2  (the  GM 
formalism)  and  Section  3  (the  similarity  argument  for  E  B(co)).  From  these  I  will  obtain 
*PV  (kz),  *PU  (cd),  4/v  (to),  Tu  (kj  and  *FV  (k^).  In  addition  it  will  be  possible  to  obtain  the 

constants  of  order  unity  such  as  a  and  a'  from  certain  consistency  requirements  and,  in 
addition,  the  quantities  kz<  and  k^will  be  obtained  from  dimensional  and  physical 

considerations. 


4. 1  The  Vertical  Wavenumber  PSD's  of  Horizontal  Velocity  ¥u  (kz)  and 
Vertical  Velocity  (kz) 

The  equation  for  *PU  (kz)  was  given  as  Eq,  (2.6).  Following  Smith  et  al  (1987)  the  value  of  a 
in  this  equation  can  be  obtained  from  the  relation 


N2  = 


k2dk  . 
z  z 


(4.1) 


This  represents  (Dewan  and  Good,  1986  and  Smith  et  al,  1987)  an  overall  saturation  condition 
(cf.  alternate  interpretation  by  Hines,  1990). 

When  the  operations  are  carried  out  is  is  possible  to  obtain  (in  a  manner  similar  to  Smith 
et  al,  1987) 
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(4.2) 


In  Dewan  and  Good  (1986)  and  Smith  et  al  (1987)  the  theoretical  value  for  k  was  not 

zmax 

discussed.  An  indication  Is  however  to  be  found  In  Dewan  (1985)  where  It  Is  argued  that  the 

largest  possible  wave-number  for  a  "wave"  Is  given  by  kmax  =  2ji  (n3 /e)  .  it  was  pointed  out 

that  this  wave  number  would  correspond  to  the  case  where  the  shortest  period  wave  would  lose 
a  substantial  amount  of  energy  through  mode  interaction  In  one  period,  thus  rendering  It  close 
to  "turbulence".  "km  "  can  also  be  recognized  as  the  well  known  "buoyancy  wave  number".  I 
have  chosen  the  2 n  following  Woods  (1974).  This  "buoyancy  wavenumber"  choice  for  kmny  has 

already  been  mentioned  In  passing  by  Munk  (1981);  but,  he  did  not  state  his  rationale.  In 
contrast,  Hines  (1990)  suggested  that  Eq.  (4.1)  be  solved  in  order  to  determine  km  . 

The  value  of  a'  In  Eq.  (3.2)  for  E  B(co)  can  be  obtained  with  the  help  of  the  GM  formalism  In 
the  following  way.  From  Sections  2  and  3  we  have 


fN 

(K)  = 

Ja>| 


A(X  )  ~2 

EB(C0)  V-S/-U  (CD)dCO 

IV 


(4.3) 


A.  m  k  /k 

z  z  z. 


(4.4) 


U  (CO)  *  (oo2  +  co2)/co2 


(4.5) 


Aog=2(i+\  f 


(4.6) 


and 


E  B(©)  -  • 

O) 


(4.7) 
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If  Eqs.  (4.4)  to  (4.7)  are  Inserted  into  Eq.  (4.3)  and  the  integral  performed  and  if  the  result  Is 
equated  to  the  expression  for  'Fu  (kz)  in  Eq.  (2.6),  then  one  can  obtain  a  relation  between  a  and 

a': 


aN2  .  / 8\  a^e 
k2  ~  V  3/  to, 

Z.  1 


(4.8) 


Solving  for  a' 


a'  =  3. 
8 


N2©. 


k2  e 

Z. 


a. 


(4.9) 


Notice  the  dimensionless  combination  In  the  square  brackets.  Let  us  set  It  equal  to  a 
dimensionless  constant  and  solve  for  k^  to  obtain 


k 


z. 


(4.10) 


where  Cz  Is  a  dimensionless  universal  constant  to  be  obtained  from  experiment.  The  physical 
meaning  behind  Eq.  (4.10)  Is  Important.  Both  Hines  (1990)  and  Munk  (1981)  have  pointed  out 
that  the  large  scale  waves  can  doppler  shift  the  small  scale  waves  and  that  for  a  certain 
critical  value  of  the  vertical  wave  number,  mc  (to  use  their  notation),  the  waves  can  be  forced 
to  approach  their  critical  layer.  Hines  gives  mc  =  N/(2  vrmg)  where  vmM  represents  the  root 
mean  square  velocity  due  to  the  large  waves.  Munk  refers  to  his  mc  as  the  boundary  for 
"compliant  waves". 

In  Eq.  (4.10),  can  be  seen  to  represent  the  quantity  mc  above  provided  that  the  "cascade" 

picture  be  kept  In  mind.  The  reason  Is  as  follows.  It  has  already  been  mentioned  in  Section  3 
that  v2  (©,)  -  e/to,  and  since  E  B(co)  ~  co-2  It  follows  that  the  main  contribution  to  v2  Is  v2  (cot) 

and  hence  It  Is  approximately  true  that  vrmg  ~  e/co^  Hence,  k^  -  N/v^  in  Eq.  (4.10)  which  Is 
to  say  that  k^  Indeed  separates  free  from  compliant  waves  as  suggested  by  Hines  (I  am  Ignoring 

constants  of  order  unity  In  this  statement). 

Next,  consider  4*v  (kz),  the  PSD  for  vertic  i  velocity.  The  GM  formalism  gives 
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(4.11) 


— 2  A(\  } 

EB(co)V  (co)  V_£/  dto 

A. 


_ 2 

where  V  (co)  = 


N2 


(4.12) 


Using  Eq.  (4.7)  for  E  B(co),  Eq.  (4.6)  for  A(XZ)  and  Eq.  (4.9)  for  a'  and  the  fact  that  N  »  colt 


¥ 


(kJ“f 


(4.13) 


2 

Comparing  Eq.  (4.13)  to  (kz)  =  a  -  (1  +  Xz)~3 

K 

z. 


(kz)  -  3  "» 

*u  (kz)  4  N 


(4.14) 


This  ratio.  It  is  easy  to  show,  also  holds  for  *FU  (k^)  and  '¥v  (k^.)  ,  and  since  a  typical  value  of 
(co,/N)  is  10-2  it  Is  clear  that  the  vertical  wave  motions  are  usually  about  two  orders  of 

magnitude  smaller  than  the  horizontal  motions.  The  value  of  a  In  Eq.  (4.13)  Is  given  by 

(1/2 
since  kmax  a  2  n  (n3 /e)  and 

K.  mCz(^  V«) ; ). 

Notice  that,  when  one  sets  cot  =  0  In  some  of  the  above  equations,  some  results  are  obtained 
which  are  not  experimentally  correct.  For  example  'Vv  (kz)  Is  presumably  not  zero  at  the 
equator.  Thus  this  model  Is  not  applicable  to  regions  close  to  the  equator. 


4,2  The  Temporal  Frequency,  co,  PSD’s:  (co),  *Py  (co) 

Using  the  GM  formalism  we  can  obtain  the  horizontal  velocity  PSD  from 
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(4.15) 


rkzm„/kz.  ^ 

E  B(co)  A(XJ  U  (co)dXz 

Jo 


~2 

Next  we  use  Eq.  12.3),  Eq.  (4.6)  (for  A(XZ)),  Eq.  (4.7)  (for  E  B(co)),  Eq.  (4.5)  (for  U  (co)),  Eq.  (4.2)  (for 
a)  and  Eqs.  (4.9)  and  (4.10)  (for  a'),  to  obtain 


*u(°>)- 


®  (#5) 


/  2  2\ 
CO  -Kl)t 


V 


CO 


(4.16) 


The  quantities  in  Eq.  (4. 16)  are  given  in  units  of  radians  (eg  radians/sec)  and  *PU  is  in 
(m2/s2)  /(radians/sec).  In  comparing  Eq.  (4.16)  with  certain  data  these  units  must  be  changed 
to  (m2/s2)  /(cydes/sec).  To  make  the  conversion  one  multiplies  the  right  side  by  2  it. 

It  is  important  to  note  that,  in  Eq.  (4.16),  when  co  »  co,,  the  frequency  dependence  goes  as 

co-2.  This  is  similar  to  the  frequency  dependence  in  the  GM  model.  The  frequency  PSD  for 
vertical  velocity  fluctuations  is  given  by 


'M®)  = 


A(Xz)EB(co)V  (co) 


(4.17) 


~2 

where  V  (co)  is  given  by  Eq.  (4.14).  In  a  manner  exactly  similar  to  Eq.  (4.16)  the  result  is 


(4.18) 


again  with  the  previous  2  n  option.  Also  the  dependence  when  co  »  tOj  is  similar  to  the  GM 
model  in  that  the  dependence  goes  as  o>°;  that  is  to  say,  it  is  constant. 


2 

t  Note  that  this  implies  that  Vy  =  ¥u  ^_when  co  »  co,. 

N2 
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4.3  The  Horizontal  Wavenumber,  kx,  PSD's:  *PU  (kx),  '?v  (kx) 

n /q  B/q 

As  has  been  mentioned  the  form  'V  (kz)  =  (CONST.)  e  kx  was  derived  in  Dewan  (1979) 
for  both  velocity  components  (u  and  v).  In  that  paper,  however,  certain  things  were  left  out  of 
account.  For  example  the  following  items  were  omitted:  (a)  the  wave  function  effects  on 
velocity  component  PSD's,  (b)  the  roll  off  effect  seen  in  the  GM  formalism  that  takes  place  at 
the  dominant  wave  number  k^,  and  (c)  the  numerical  value  of  the  constant  of  order  unity.  The 

purpose  of  the  present  section  is  to  finish  this  task. 

From  the  above  mentioned  form  for  'P  (kj  one  is  led  to  the  assumption,  in  the  context  of 

the  GM  formalism,  that 


where  Xx  *  k^/k^.  I  have  chosen  to  commence  with  the  horizontal  velocity  component  because 
the  latter  is  involved  in  the  Hines-Munk  critical  layer  mechanism.  From  the  GM  formalism 


E  B(<a) 


~2 


U  (o>)dco. 


(4.20) 


_fi/q 

Again,  A(X)  =  (s  -  1)(1  +  X.)-"  and  setting  s  =  5/3,  A  (XJ  =  (3)  0  +  K)  •  AS®1”  using 
E  B(io)  «  a'  e/a >2,  and  U  (<o)  =  (co2  +  ©2)/<o2  and  performing  the  integral, 


k  o), 

x.  1 


(4.21) 


The  value  of  a"  is  obtained  by  equating  Eq.  (4.21)  to  Eq.  (4.19) 


£l/3k2/3N 

- ^  /r)  a' 

cot  J  \9/ 


(4.22) 
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from  Eqs.  (4.9)  and  (4.10)  and  Eq.  (4.23)  below,  t 

Again  we  note  the  dimensionless  combination  in  Eq.  (4.22)  and  in  a  manner  similar  to 
Eq.  (4.10)  we  obtain 


r~r 


(4.23) 


which  has  an  important  physical  interpretation  in  the  context  of  the  Hines-Munk  critical 
layer  picture.  The  horizontal  phase  trace  velocity  is  given  by  co/kx.  For  a  critical  layer  to 

occur.  (to/kx)  -  urma  where  k^  Is  the  critical  wave  number  at  that  frequency,  ©.  Unlike  the 

case  for  kz  where  only  one  kz  Is  critical,  one  has  k^  (co)  which  Is  to  say  It  is  frequency 

dependent.  From  the  cascade  assumption  (e/coj  -  u^ms  hence 


„  2 
_2  <fl  CD 

K  — r1  (4.24) 


Thus  the  physical  meaning  of  k^  is  that  it  is  the  value  of  kx  for  the  lowest  wave  frequency.i.e. 
to  =  ©,. 

It  should  be  noted  that  since  kx  (©)  depends  on  ©,  this  may  provide  a  clue  to  the 

Interpretation  of  the  ©  dependence  found  by  Pinkel  (1981,  1984).  That  is  to  say,  when  a 
*u  (®*  kx)  spectrum  Is  obtained,  our  approach  would  indicate  that  at  frequency  ©, 

and  that  the  dominant  wave  number  is  frequency  dependent.  That  is  to  say 
the  roll  off  point  in  k^  dependence  would  change  with  frequency. 

t  As  has  been  mentioned,  when  Eq.  (3.2)  or  (4.7)  is  integrated  between  ©,  and  N,  and  the 

normalization  of  B  (©),  see  (2.3),  is  used,  one  obtains  E  =  a  e  ©~  .  But  E  =  q2,  the  velocity 
variance;  thus,  to  within  a  constant  of  order  unity,  a',  we  have  q2  =  e/©r 
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In  order  to  obtain  '¥v  (1^)  one  makes  use  of  Eq.  (4.14)  which  also  holds  for  *FU  v*  (k^).  That 
Is  we  use  '¥y  (kx)  =  |^-Tu  (kj  obtaining 


4» 


(kx)  = 


*5 

4  N 


,2/3 


e2/3 
i  5/3 


(1+\) 


-5/3 


(4.25) 


A  final  caution  Is  In  order  here.  In  the  case  mentioned  In  Dewan  (1979)  the  fluctuations 
were  exceptionally  severe  and  were,  In  fact,  previously  Interpreted  as  ordinary  turbulence.  In 
addition,  the  velocity  components  were  comparable  In  size.  The  present  model,  however,  Is 
not  Intended  for  unusually  high  energy  situations.  It  Is  Intended  to  represent  the  usual 
"background"  "universal”  case.  In  this  sense,  the  experimental  data  discussed  In  Dewan  (1979) 
do  not  represent  "confirmation"  of  the  present  model.  The  latter  awaits  future  experiments. 

A  similar  case  of  exceptionally  high  energy  gravity  waves  Is  discussed  In  VanZandt  et  al 
(1990)  where  *PV  (co)  falls,  not  as  <o°,  as  In  the  less  excited  case,  but  as  (as  would  be 
expected  for  4*u  (co)).  In  both  of  these  situations,  (that  Is,  Dewan  (1979)  and  VanZandt  (1990)1 
the  observed  departure  from  the  predicted  PSD  may  be  due  to  the  cause  suggested  by  Van  Zandt 
(1990).  This  is,  namely,  that  the  vertical  velocity  component  measurements  are  contaminated 
by  the  effects  of  the  very  large  horizontal  components  of  the  velocities. 


5.  COLLECTION  OF  PREDICTIONS 


5.1  Model  Equations 


^u(kz)  =  ( 1  +  K)~3  (m2/s2)/  (rad/m) 

k 

z. 

x,*(Vk0-  vc.VK®,/') 

“-[ln(kv„/k0- (I)]"' 


1  This  means  Tu  and  *PV. 
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(5.2) 


¥ 

V 


"Mr"  ( 1  +  \f  (m2/s2)/(rad/m) 
k 

Z. 


2/3 


^uiK)  =^rh3~(1+K)'*/3  (m2/s2)/(rad/m) 


k 

aC 


5/3 

,2/3 


a  = 


3C 


2  ’ 


\  *k  /k 

X  X  X 


k  =cxV(cof/e) 


(5.3) 


(5.4) 


^(co)  = 


a'e 

2 

(O 


2  2  "I 

co  +  a>j 


(2n)(m2/s2)/  (cydes/sec) 


(5.5) 


*v(«) 


a'e 


2 

co 


N2 


(2  it)  (m2/s2)/  (cydes/sec) 


(5.6) 


5.2  Suggestions  for  Future  Tests  of  This  Model 


Clearly,  future  tests  would  Involve  measurements  of  e  simultaneously  with  power  spectral 
measurements  of  velodty  fluctuations.  In  situ  and  remote  measurements  of  e  have  been 
performed  by  Barat  (1975),  Cadet  (1977),  and  Hocking  (1988a,  1988b,  1985)  for  example. 

The  values  of  k^  and  k^  could  be  obtained  directly  from  the  relevant  PSD's  (from  the  "roll 

off'  characteristics)  and,  by  us»ng  known  values  of  N,  e,  and  co,,  the  values  of  Cx  and  Cz  can  be 
obtained  from  the  equations  for  kx>  and  kz  above.  The  dependence  of  these  and  the  PSD's  on 

co,  could  be  checked  by  conducting  experiments  at  various  latitudes.  Such  dependences  are  all. 
of  course,  entirely  new  predictions. 
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It  seems  that  the  most  convenient  tests  of  the  predicted  e  dependence  of  these  gravity  waves 
would  Involve  the  frequency  PSD's  (Eqs.  (5.5)  and  (5.6)).  If  these  obtain  verification  It  would 
then  be  appropriate  to  test  the  more  difficult  to  measure  'FfkJ  PSD's  for  e  dependence. 

Oceanic  measurements  of  e  together  with  appropriate  PSD’s  would  also  play  an  Important 
role  In  the  verification  of  this  model  In  spite  of  the  difference  of  N  dependence  upon  depth. 
Sonar  would  then  replace  radar  as  the  Instrument  of  choice. 


6.  ORDER  OF  MAGNITUDE  COMPARISON  WITH  EXPERIMENT 


If  simultaneous  measurements  of  PSD's  and  e  in  the  atmosphere  or  ocean  were  presently 
available  it  would  be  possible  to  put  the  present  model  to  a  test  immediately.  Unfortunately, 
no  such  data  seem  to  exist  yet.  It  Is  nevertheless  possible  to  find  out  If  the  theory  predicts 
quantities  "within  the  bounds"  of  observation.  In  carrying  out  this  exercise  I  learned  that  the 
model  Is  rather  elaborately  "Interconnected"  and  hence  "vulnerable"  to  quantitative 
observation.  For  example,  k2  and  Cz  enter  Into  the  PSD  coefficient  In  both  a  direct  and 

(through  a)  Indirect  way. 

In  the  following,  "reasonable"  values  for  the  unknowns  were  used  and  the  relevant  papers 
cited. 


6. 1  Atmospheric  Comparison 

As  mentioned  in  Section  4.1,  the  *FU  (kz)  PSD  for  the  atmosphere  has  been  found  to  be 
valid.  We  therefore  turn  to  a  frequency  PSD  namely  ¥  (to)  and  compare  It  with  experiment. 
For  convenience  I  list,  on  the  basis  of  Section  (5.1),  the  following  useful  equations: 


k 


X  »2x/k  , 

*.  Z. 


X 

zroaj 


2n/k 

zmax 


(6.1) 


(6.2) 


VJN)  =  (|)  C'2£  (l/N*)2  na 


(6.3) 
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“-iln(k.„„/k0- (§)]"' 


(6.4) 


The  "typical"  values  for  the  parameters  will  be  chosen  to  be 


c  = 

5  x  10”4  m2/s3 

Cadet  (1977) 

N  = 

1CT2  rad/s 

VanZandt  (1982),  his  Fig.  1 

to,  = 

10-4  rad/s  (~40°N) 

VanZandt  (1982),  his  Fig.  1 

Cz  " 

0.2 

(chosen  because  of  experimental  agreement) 

a  = 

0.235 

[from  (6.4)J 

From  these  the  following  results  can  be  obtained:  k  ~  8.9  x  10^*  rad/m,  X_  =  7  km, 

Z*  Z* 

k.  s  2.8  x  10-1  rad/m,  X„  =  22  m.  A  comparison  between  these  numbers  and  those  found 

zmax  zmln 

in  Smith  et  al  (1987)  (X  in  stratosphere  «>  5  km),  and  Dewan  and  Good  (1986)  (X  ~  40  m) 

zmln 

shows  that  these  length  scales  are  acceptable.  The  PSD's  Tu  (to,)  and  *P  (N)  are 
S*u(®i)  =  1-4  x  106  (m2/s2) /(cycles/sec);  [cf.  105  -  106  3] 


Tu(N)  =  69  (m2/s2) /(cycles/sec);  [cf.  101  8  or  63  -  103  5  . 


The  values  indicated  by  "cf."  are  the  experimental  ranges  taken  from  VanZandt  (1982), 

Fig.  1.  Thus,  I  conclude  that  if  I  choose  Cz  =  0.2  I  can  get  "good"  agreement  with  experiment. 

Is  the  value  for  Cz,  "a  universal  constant  of  order  unity",  a  reasonable  value?  The  answer  is 

yes  as  can  be  shown  by  considering  Bond  (1929). 

It  has  been  mentioned  that  VanZandt  (1982)  observed  that  Tu  (to)  had  a  to~5/3  rather  than  a 

to-2  dependence.  In  this  case  one  would  be  forced  to  alter  the  present  model  along  the  lines  of 
Tcnnekes  (1975)  as  was  mentioned  earlier  in  Section  3.  See  Eq.  (3.3).  This  modification  is 

very  simple.  E  B(co)  would  then  become  proportional*  to  f — ^773!  ®  .  I  feel  that  such  a 

l(®i)  ) 

modification  would  be  indicated  once  the  present  model  has  been  at  least  partially  validated 
by  direct  measure-ments  of  ewhen  XFU  (to)  or  *FV  (cn)  were  simultaneously  measured,  and  if  the 

of5/3  continued  to  be  observed. 


*  A  detailed  argument  for  why  to  5/3  should  be  correct  for  the  atmosphere  while  to-2 
applies  to  the  ocean  will  be  found  in  the  appendix. 
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It  is  possible  to  put  the  present  model  to  one  further  test  at  this  time.  According  to  Eq.  (5.5) 
the  (co)  PSD  is  linearly  dependent  on  e.  Measurements  of  e  in  the  mesosphere  have  been 

made  by  Hocking  (1988)  and  found  to  be  of  order  0.  lm2s-3  at  86-88  km  altitude.  Comparing 
this  to  e  =  5  x  10-4m2s-3  found  in  the  stratosphere  by  Cadet  (1977)  (1 1  km)  it  is  found  that  the 
present  theory  predicts  that  Tu  (co)  at  the  higher  altitude  can  be  200  times  as  large  as  it  is  at 

the  lower  altitude.  Has  this  been  observed?  According  to  Balsley  and  Carter  (1982)  "In  the 
co"5/3  portions  of  the  spectra,  the  spectral  energy  density /unit  mass  is  about  250  times  smaller 
in  the  troposphere  (8  km)  than  in  the  mesosphere  (86  km)".  This  is  encouraging.  Of  course  a 
valid  test  would  require  simultaneous  measurements  of  e  and  *FU  (co)  at  such  altitudes.  For 
example,  Balsley  and  Garello  (1985)  found  a  ratio  of  about  50  in  'P  (co)  between  13  km  and 
87  km. 

Between  50  and  250  we  have  a  factor  of  5;  but  the  GM  model  claim  is  that  these  PSD’s  are 
generally  within  a  factor  of  2.  This  might  imply  that,  in  this  sense,  the  atmosphere  is  more 
variable.  An  interesting  discussion  about  e  and  universality  will  be  found  in  Munk  (1981), 
sections  9.9.6  and  9.9.7. 

It  can  be  seen  from  the  above  that  in  the  case  of  atmospheric  spectra  such  as  *P  (co)  and 
*u  v^Sc  the  present  model’s  dependence  on  e  makes  it  desirable  to  abandon  the  use  of  the 
word  "universal"  for  such  PSD’s.  Since  e  varies  over  a  large  range  in  the  ocean  it  is  curious 
that  these  PSD’s  have  not  been  reported  as  being  very  variable  there. 


6.2  Oceanic  Comparison 

The  appropriate  values  for  "typical"  oceanic  values  are  as  follows: 


e  * 

KT8  (nr2  /sec3) 

Woods  (1974)  and  GM  (1972),  p.  251 
also  Monin  and  Ozmidov  (1985) 

/s 

N  = 

3.5  x  KT3  (rad/ sec) 

GM  (1972),  2  cycles/hv 

(0  = 

1 

7  x  10“5  (rad/sec) 

GM  (1972),  0.04  cycles/hv,  p.  261 

M  = 

7.7  x  KT4  (rad/m) 

GM  (1972),  0.122  cycles/km 

E'  = 

2  jix  10-5 

GM  (1972),  p.  252 

Cz  - 

(0.2) 

(from  previous  section) 

Next  to  each  quantity  at  least  one  reference  is  given  for  support.  Perhaps  the  most  variable 
parameter  Is  e  which  Monin  and  Ozmidov  (1985)  indicate  can  vary  over  a  very  wide  range 
depending  upon  depth,  location,  and  weather  conditions.  Values  listed  by  them  (pp.  124-129) 
go  from  1(T10  to  10"5  m2/sec3.  The  "hat"  notation  is,  as  was  mentioned  above,  from  GM  who 
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use  it  to  mean  "dimensional"  to  distinguish  the  quantities  in  their  model  which  are  primarily 
non-dimensional.  The  quantity  Cz  was  set  equal  to  0.2  as  it  was  in  the  previous  atmospheric 
comparison  since  Cz  is  presumed  to  be  a  universal  dimensionless  constant. 

As  has  already  been  mentioned  the  GM  model  gives  a  scaling  with  respect  to  N  which 
differs  from  the  present  one.  This  matter  cannot  be  further  investigated  here.  In  what  follows 
one  can  only  determine  if  the  present  model  is,  in  some  sense,  reasonable.  The  following  are 
presented  for  comparison: 


GM  model: 


r 

2  ^1 

r  -2  ^2  i 

|  2E' 

f  N  ^ 

C0  +  COj 

L  * 

Im; 

J 

~3  J  ^2  ^2 

(0  V  CO  -  C0j 

(6.5) 


Present  model:  4^(00)  =  _^(2n)  ^3^ 


r  2  2 

(0  +  COj 


0) 


(6.6) 


o*[ln(kw/k0‘‘-T' 


and 


GM  model:  4*  (k^)  » 


^(1.5)  (3/2) 

k 

Z. 


-2.5 


(6.7) 


Present  model:  4*  (kz)  = 


(6.8) 


The  values  k,  ,  ,  k  „  and  are  obtained  from  Eq.  (6.1).  Table  6.1  shows  the 

comparisons  between  the  models. 


Table  6. 1  Comparison  of  Present  Model  with  "Observations". 


Present  Model 

X„  =  0.48  m 

zmln 

X,  =  107  m 

z» 

a  =0.256 

=  5.9  x  10-2  rad/m-1 

k  =  13  nT1 

zmax 

"Observation" 

X  =~  1  meter  Woods  (1974) 

Tnln 

X^=  170  m  Holloway  (1981),  p.  58 

(2rt)  (|)  ~  =  1.51  x  10'7  (m2/s2)  (E|d) 

L  Z  J 

=  part  of T  (co)  [Eq.  (6.6)] 

^~(§)  =5-8xl0"6(m2/s2)(E|Sl) 

=  part  of  ^  (o)^  [Eq.  (6.5)] 

aN2  ,  „  ._-2  m2/s2 

— —=1.6x10  - — -i- — - 

k3  (rad/m) 

z. 

=  partof'Pu(kz)  [Eq.  (6.8)] 

"  E'  =  5 x  10~2 
~2  kz.  (rad/m) 

=  partof'Fu(kz)GM 

_ [Eq-  (6-7)] _ 

f  In  order  to  avoid  problems  with  dimensional  inconsistency,  it  would  appear  necessary 

/  2 

to  ascribe  the  units  — y  to  e  in  the  present  context  (pointed  out  to  me  by  T.  Van  Zandt 

V  ®  / 

in  a  private  communication). 

*P  (©j)  =  62  (m2/s2) /(cycles/sec) 

[Eq.  (6.6)] 

56  (m2/ s2) /(cycles/sec) 

*P  (N)  =  1.2  x  10~2  (m2/s2)/(cycles/sec) 

[Eq.  (6.6)1 

1.7  x  10-2  (m2/s2)/(cycles/sec) 

Both  from  GM  (1972),  Fig.  4.  See  text. 

Note  that  the  4,u  (<d)  comparisons  were  performed  in  three  ways.  First  the  values  of  the 
"curly  brackets"  in  Eqs.  (6.5)  and  (6.6)  were  compared,  then  *PU  (co,)  from  Eq.  (6.6)  was  compared 
with  the  value  derived  from  Fig.  4  of  GM  (1972).  The  same  was  done  for  4*u  (N)  from  Eq.  (6.6). 
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To  make  the  comparisons,  the  graphical  values  had  to  be  multiplied  by 

3.5  x  lCT3  rad/  s  r;  .  ,  . 

- - —  =  N  In  cycles/sec. 

2xrad/cy 

Regarding  (k2)  there  was  a  difference  In  slope.  The  present  model  gives  (1  +  X.z)-3  In 
comparison  to  (1  +  Xz)-2-5  InGM  (1975).  We  do  not  know  whether  or  not  the  difference  is  real  or 


due  to  a  difference  in  data  processing.  For  example,  omission  of  "pre-whitening"  techniques 
could  explain  the  differences.  As  can  be  seen,  the  values  of  the  curly  brackets  are  comparable, 
In  general,  the  numerical  values  in  Table  6. 1  agree  to  within  a  factor  of  about  3  or  better. 

What  has  been  proven  in  Sections  (6.1)  and  (6.2)?  All  that  was  accomplished  was  that  the 
model  appears  to  give  the  correct  order  of  magnitude.  In  some  circles  this  Is  called  a  "sanity 
test".  The  model  is  not"obviously  wrong",  In  other  words. 


7.  CONCLUSIONS 


A  similitude  model,  based  mainly  on  a  wave  interaction  cascade,  has  been  given  to 
"explain"  the  GM  model  from  physical  considerations.  By  necessity,  the  model  has 
experimental  consequences  that  remain  to  be  tested.  These,  at  present,  are  Its  main  virtue.  If 
future  experiments  validate  the  model  then  it  will  have  demonstrated  predictive  value.  The 
main  experiments  suggested  include  simultaneous  measurements  of  e  and  certain  PSD's  as 
given  in  Section  6. 

The  exact  roles  of  "saturation"  due  to  "Instability"  and  "cascade"  due  to  "nonlinear  mode- 
interaction"  are  not  entirely  clear  at  present;  but  it  appears  as  If  both  are  taking  place. 
Holloway  (1980)  has  pointed  out  that  nonlinear  mode- Interaction  can.  In  and  of  Itself, 
represent  "saturation"  of  a  type  with  critical  Richardson  number  «*  4.  Both  types  of 
"saturation"  might  take  place  and  be  related  to  each  other  and  bring  about  the  observed 
spectra.  Obviously  the  situation  is  not  completely  understood  yet  however. 

As  has  been  mentioned  previously  there  has  been  some  question  of  whether  or  not  the 
waves  discussed  in  this  paper  should  be  called  turbulence  of  some  kind  (e.g.  "buoyancy 
turbulence",  or  "2-dimensional  turbulence")  -  see  Gage  (1979),  Trefethen  and  Panton  (1989),  Sidi 
et  al  (1988),  Weinstock  (1985),  Holloway  (1983)  and  Gargett  et  al  (1981).  The  above  discussion 
shows  that  wave  properties  exist  (like  propogation,  dispersion  relations,  wave  functions,  and 
polarization  relations)  and  yet,  at  the  same  time  there  are  some  of  the  properties  of  turbulence 
such  as  self-similar  cascades  In  these  velocity  fluctuations.  We  seem  to  be  dealing  with  a 
hybrid  dynamics  for  which  no  name  yet  exists.  Does  this  mean  that  we  need  a  new  word  like 
"waveulence"  or  "turbundulence",  or  will  "buoyancy  turbulence"  be  altered  in  meaning  to  stand 
for  the  above? 

Regarding  the  term  "universal  spectrum",  the  predicted  dependence  of  4*u  v  (to)  and  4*  v  (kx) 
upon  e  Indicates  that,  If  the  model  proposed  here  is  valid,  these  spectra  would  be  much  "less" 
universal  than  the  ¥u  v  (kz)  spectra.  It  was  also  mentioned  that  the  slope  of  the  VPU  (o>) 
spectrum  may  be  -2  or  -5/3  depending  upon  the  conditions  (see  the  Appendix). 
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Appendix 


Convective  Effects  upon  the  Eulerian  Spectrum  for  Energy  as  a  Function  of  to:  EB(co) 


In  the  text  it  was  menUoned  that  the  frequency  spectrum  of  horizontal  velocity 
fluctuaUons  had  an  co~2  dependence  for  internal  gravity  waves  in  the  ocean  whereas  in  the 
atmosphere  van  Zandt  showed  evidence  for  an  uf5^3  dependence.  It  was  also  mentioned  in  the 
text  that  Tennekes  (1975)  gave  an  argument  in  the  context  of  turbulence  that  might  possibly  be 
adapted  to  the  case  of  waves  and  might  therefore  be  useful  in  explaining  the  difference  in  oo 
dependence  in  oceanic  and  atmospheric  observations.  The  purpose  of  this  appendix  is  to  carry 
out  this  suggestion. 

Tennekes  (1975)  pointed  out  that  the  Kolmogorov  dimensional  arguments  for  inertial 
range,  or  cascade,  turbulence  lead  to 


(A- 1 ) 


where  <|>L  (co)  is  the  velocity  fluctuaUon  PSD,  the  subscript  L  designates  "Lagrangian 
descripUon",  and  PL  is  a  constant  of  order  unity.  He  then  pointed  out  that  experimental 
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evidence  contradicts  (A- 1 ) .  He  next  showed  that  for  sufficiently  large  Reynolds  number  based 
on  "outer  length",  the  inertial  range  would  turn  Into  an  "Inertial  advective"  range  which 
Includes  the  advective  effects  of  large  eddies  upon  the  small  eddies.  This  changes  the  Eulerian 
spectrum  which  now  would  no  longer  resemble  the  Lagrangian  spectrum  due  to  the  doppler 
effects.  More  specifically  he  showed  that  advective  doppler  effects  cause  the  Eulerian 
frequency  PSD  for  velocity  fluctuations  to  have  the  following  form: 


0.  (o>)  *  P, 


(q2/3£2/3) 


0) 


5/3 


(A-2) 


where  Pg  is  a  constant  of  order  unity,  the  subscript  E  designates  Eulerian  frame,  and  q  is  the 
root  mean  square  velocity  Auction  due  to  the  turbulent  field.  This  relation  will  be  derived 
below. 

In  order  to  ascertain  the  conditions  necessary  for  advective  effects  to  become  important, 
Tennekes  made  use  of  the  fact  that  such  effects  are  most  prominent  at  the  highest  frequencies. 
In  inertial  range  turbulence  the  highest  frequency  is  known  to  be 


0) 


c.o. 


(A- 3) 


where  the  subscript  c.o.  refers  to  "cut  off'.  The  kinetic  energy  {per  unit  mass)  in  the 

1  2  1  2 

Lagrangian  and  Eulerian  frames  at  coc  Q  are  designated  ^-u  L  (coc  o  )  and  ~  u  g  (g>c  q  ) 
respectively.  Note  that  coc  0  is  the  highest  frequency  for  the  L  frame  but  not  for  the  E  frame 
since  doppler  effects  can  make  frequencies  shift  upwards.  Tennekes  compared  these  energies 
and  found 


2  U  E  (“c.o.)  .  1/6 

* - =  Re, 

J-uf(o)  ) 

2  LV  c.o./ 


(A-4) 


where  1  is  the  outer  scale  of  the  cascade  and  Re  is  the  Reynolds  number.  Thus  for  large  enough 
Re,  the  Eulerian  energy  at  the  Lagrangian  frequency  cut  off  can  be  significantly  larger  than  the 
corresponding  energy  in  the  Lagrangian  frame  due  to  advection. 

In  order  to  deduce  the  PSD  in  the  Eulerian  (E)  frame,  that  is  (A-2),  Tennekes  used 
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(A-5) 


which  Is  the  frequency  in  E  due  to  q  advective  doppler  effects  on  eddies  of  scale  r  caused  by  the 
largest  scale  eddies.  From  Kolmogorov's  theory  he  used 


o  ul(°>)  ~  e 


.2/3  r2/3 


(A-6) 


Inserting  (A-5)  into  (A-6).  by  elimination  of  r. 


lu  !(<*>)  ~e2/3 


(A-7) 


Using  the  fact  that  the  PSD  is  the  "kinetic  energy",  i.e.  ^u 

£ 


per  unit  frequency  he  obtained 


_  PE  (eq) 
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CO 


5/3 


QED. 


(A-8) 


This  differs  from  (A- 1)  in  the  same  way  that  atmospheric  waves  seem  to  differ  from  oceanic 
waves.  Thus  in  the  case  of  cascade  turbulence  (A-2)  is  valid  when  Re  in  (A-4)  is  sufficiently 
large. 

Now  we  turn  to  the  question  of  what  happens  to  the  frequency  spectra  of  waves.  As  will  be 
shown,  a  "frozen  waves"  hypothesis  which  is  analogous  to  the  "frozen  turbulence"  implied  by 
(A-5)  appears  to  be  justified.  I  shall  start  by  assuming  it  is  valid  and  then  subsequently 
present  evidence  that  it  is  indeed  justified. 

We  start  by  deriving  a  counterpart  to  (A-4)  for  waves  with  the  idea  that  the  result  will 
possibly  indicate  a  difference  between  ocean  and  atmosphere.  The  highest  or  "cut  off' 
frequency  for  internal  waves  is  the  buoyancy  frequency  N.  This  replaces  the  oc  0  in  the 

Tennekes  argument.  The  first  question  to  answer,  then,  is  what  is  the  value  of  ^  u  £  (N)? 

From  the  main  text  we  have 


4>L(w)aEB(co)  =  ^ 
(0 


(A-9) 
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Setting  co  =  N  In  (A-8)  and  Aco  =  co  as  discussed  In  the  text  In  the  context  of  fractal  self- 
similarity,  we  obtain 


2  L  u  N 


(A- 10) 


In  order  to  obtain  i  u^(N)  we  follow  Tennekes  by  assuming  the  frozen  hypothesis  and 
z 

commence  with  Lagrangian  spatial  dependence.  From  the  main  text  we  have 


2/3  .  -5/3 


Vu  (kz)  =  e  k 


(A-l  1) 


Notice  here  that  I  have  chosen  the  horizontal  velocity  components  only.  This  is  because 

the  vertical  components  at  any  wave  number  are  smaller  by  a  factor  of  j  as  shown  in 
Equ.'s  (4.14)  and  (5.4).  Next,  the  frozen  assumption  gives 


=  <1  k 


(A- 12) 


which  is  our  counterpart  to  Tenneke's  (A-5)  above.  This  assumes  that  for  the  high  Eulerian 
frequencies,  coE,  one  senses  only  the  convective  effects  due  to  spatial  scales  which  are  caused 

by  the  wave  field  as  a  whole.  The  latter  effect  is  of  course  due  to  q,  the  r.m.s.  velocity  (compare 
Hines,  1991).  Solving  (A-12)  for  k^  and  substituting  the  result  Into  (A-ll)  we  obtain  (under  our 

assumption)  an  estimate  for  4>E  (wE),  the  Eulerian  PSD  fcr  velocity  fluctuations: 


vuK) 


(A- 13) 


In  order  to  obtain 


(N)  we  again  multiply  by  the  appropriate  band  width. 


But  since 
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(Eq.  A- 13)  was  obtained  from  (A- 11)  which  Is  a  density,  not  In  co  but  in  k,  we  must  multiply  yu 
by  Ak.  Also,  coE  =  N  is  the  relevant  frequency.  What  then  does  one  use  for  Ak?  The  answer  is 

Ak  =  (using  A- 12  with  coE  =  N  and  using  Ak  =  k  as  before).  Thus  we  have 


A  u  2  (N)  =  »(/u  (N)  •  Ak  =  a"  e2/3 
2  ^ 


(K 

u 


n 


-5/3 


=  a  e 


2/3 


(A- 14) 


We  finally  arrive  at  the  ratio  [in  analogy  with  (A-4)l 


2^1iN)_  ,  «"  (N/q)~2/3  =  (<C)  N'/3q2/3 

iuf(N)  oc'eN'1  Va'J  t1'3 

2  L 


(A- 15) 


where  (A- 10)  was  used.  But,  as  shown  in  the  main  text. 


rN 

q2  =  I  EB  (to)  d((D)  = 
J  COj 


ae 

co 

i 


(A- 16) 


Solving  (A- 16)  for  e  and  inserting  into  (A- 15), 


1u3(N) 


(A- 17) 


to  within  constants  of  order  unity. 

This  equation  states  that  as  the  ratio  N/co(  increases,  the  convective  effects  upon  the  higher 
frequencies  increases  and,  as  a  consequence  of  the  earlier  arguments  of  Tennekes  that  were 
described  above  and  which  continue  to  hold  under  the  present  assumptions,  the  PSD  for 
energy  changes  from 
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(A- 18) 


EB  (CD,  )  ~  i 
CO 


to 


,  ,  2/3  2/3  -5/3 

EB  (coE)  ~  e  q  '  co 


_  e 


-5/3 


CD. 


1/3 


CO 


(A- 19) 


where  I  have  omitted  constants  of  order  unity.  In  Eq.  (A- 19)  I  used  (A- 14)  with  coE  replacing  N. 

The  Eulerian  PSD  was  then  obtained  by  then  dividing  ~  u  E  (o>E)  by  the  bandwidth  co£  In  the 

usual  way.  and  thus  I  obtained  the  middle  term  of  (A- 19). 

The  next  question  is  whether  or  not  (A- 17)  can  explain  the  difference  between  the  slopes 
observed  between  oceanic  and  atmospheric  temporal  frequency  spectra.  In  the  main  text, 
when  comparisons  were  made  between  observations  and  theory,  approximately  the  following 
values  were  employed. 


Atmospheric:  N  =  2.2  x  10  2  rad/sec*,  co,  =  1 0-4  rad/sect  (50°N) 

Oceanic:  N  =  3.5  x  10~3  rad/sec,  co,  =  7  x  10“ 5  rad/s  (29°M) 

(9  x  10'5  rad/s  (39°N) 


A  closer  look  at  Table  1  of  GM  (1972)  p.  237  shows  that  his  Fig.  1  is  really  based  on  data 
such  that  39° N  is  a  more  appropriate  latitude  [this  is  based  on  a  weighted  average  of  co,]  thus 
we  will  replace  co,  above,  for  the  ocean,  by  9  x  10"5  rad/s  (in  brackets).  Tills  would  not  change 
any  conclusions  in  the  main  text.  Thus,  from  (A- 16)  we  have 


t  These  values  are  not  exactly  the  same  as  those  in  Section  6.1  of  the  main  text.  There  I  used 
N  =  10-2  rad/s  and  co,  =  10-4  rad/s  (40°N).  The  present  values  here  I  believe  to  be  more 
accurate.  They  change  no  conclusions  of  the  main  text  however.  The  new  value  of  \  =  3.5  km 
instead  of  7  km,  ^.m,n  =  7.9  m  instead  of  22  m  and  \yu(N)  =  14  (m2/s2)/(c/s)  instead  of 

69(m2/s2)/(c/s)  .  This  last  value  takes  it  lower  than  the  experimental  minimum  of  63;  but, 
remember  that  the  point  of  the  present  appendix  is  to  show  that  the  -5/3  slope  Is  valid  and 
this  would  change  yu(N)  to  a  higher  value.  For  these  reasons  the  conclusions  of  the  main  text 
were  left  unchanged. 
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Atmospheric: 


u. 


6.0 


Oceanic: 


3.4 


These  are  in  the  ratio  of  (6/3.4)  =  1.8  which,  as  can  easily  be  shown,  is  just  sufficient  to 
convert  the  atmospheric  PSD  slope  from  <o~5/3  to  o>-2.  In  other  words,  reducing  the  PSD  at 
a)  =  N  in  van  Zandt's  Figure  1  by  80%  steepens  his  slope  from  -5/3  to  -2  on  the  log-log  graph. 

In  conclusion  it  appears  that  we  have  a  quantitative  explanation  of  the  difference  between 
atmospheric  and  oceanic  measurements  of  yu  (to).  One  interesting  consequence  of  this  is  the 
following  prediction.  When  (N/cOj) 1  =  6.0  in  the  ocean,  then  one  should  also  see  an  aT5/3 

dependence!  What  latitude  gives  the  appropriate  <0j?  Using  N  =  3.5  x  10“  *  rad/s,  we  are  led  to 

<0i  =  N  =  i  6  x  10~5  rad/s  (a-20) 


But  since 


co(  =  2£2sin  <)>  =  (1.46  x  1 0~ 4)  sin  d> 

(l.6  x  10-5) 

( 1.46  x  10-4) 


,  .  (1.6X10"5)  _„oNI 

<>  =  arc  sin  ------  ;  =  6.3  N 


(A-2 1 ) 


This  latitude,  unfortunately  is  very  close  to  the  equator.  As  was  mentioned  in  the  text,  the 
present  model  is  not  applicable  at  <J>  =  0°.  This  therefore  raises  the  question  of  whether  or  not 
It  Is  valid  to  predict  the  above  result  from  the  model. 

An  alternative  test  would  be  to  see  if  there  were  an  extreme  northern  latitude  where  the 

1/3 

atmospheric  PSD  would  change  to  an  oT2  dependence.  That  would  require  f  N  )  =  3.4  and 

for  N  =  2.2  x  10  2,  cof  would  have  to  be  5.6  x  10"4  rad/s  which  is  not  possible.  Perhaps  with 


35 


sufficiently  accurate  PSD's  the  implied  changes  in  PSD  slope  with  latitude  could  be  detected 
but  such  accuracy  may  be  too  expensive  to  be  immediately  practical.  The  other  tests  in  the 
text  are  much  more  practical  to  perform.  If  these  tests  turn  out  to  be  positive,  then  possibly,  a 
latitude  test  may  be  tempting  to  perform. 

We  now  turn  to  the  question  of  the  validity  of  the  "frozen  wave  assumption"  made  above. 
The  relation  between  the  Eulerian  and  Lagrangian  frequencies  where  the  Lagrangian  frame 
has  a  velocity  q  (presumed  horizontal  for  previously  given  reasons)  is 


a>E  =  <oL±qkx 


(A-22) 


2 

In  order  to  ascertain  the  value  of  uE  (N)  we  may  use 


^u(kx-a,L)d<°L*Akx 


(A-23) 
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(A-24) 


HK) 


which  represents  EB(co)  U  (co)  -V  ;c-  from  the  main  text.  When  (A-22)  is  used  to  eliminate 

k 

X. 


from  (A-24)  one  can  obtain,  lor  the  case  of  co-.  =  N  and  hence  k  =  f - ]  from  (A-22), 

x  V  8 


ue(n)=[  q3^  y°‘-)  (2/3) 

Ju),  0), 


t  C0(  i 
V  CO,  +  N  -  coL  f 
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dco.  •  Ak 
L  x 


(A-25) 


where  the  approximation 
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(A-26) 


co. 


was  used.  Notice  that  the  relative  sign  between  kx  and  q  was  chosen  as  positive,  and  therefore, 
a  positive  change  of  coL  in  the  integration  of  (A-24)  will  result  in  larger  rather  than  smaller 
contributions  from  what  was  the  kx  term  in  (A-23).  We  regard  the  contributions  of  the 
opposite  case  to  be  comparatively  negligible.  Now  Eq.  (A-24)  can  be  written 


u|(N)  =  Akx  q3  cof/3-  (INY) 


(A-27) 


where 


(I  NT)  s 


f 


M-?) 


°>L  (01,  +  N  -  (Oj 


5/3 


1  dco. 


(A-28) 


where  the  (2/3)  factor  is  now  dropped  since  Cz,  and  a'  have  already  been  set  equal  to  unity. 
Using  the  usual  order  of  magnitude  values,  co,  =  10  and  N  =  10  ~2  (rad/sec)  a  numerical 

Integration  of  INT  is  found  to  give  (3.7  +  .4)  x  107. 

But  we  seek  a  convenient  estimate  of  (A-27)  which  gives  the  dependence  on  co,  and  N 

explicitly.  It  is  clear,  by  inspection,  that  the  greatest  contribution  of  the  integrand  would  be 
for  values  of  coL  near  cor  For  example,  if  the  integrand  of  (A-27)  is  caked  f  (ooL)  [  letting  co,  and 
N  assume  the  above  values)  and  plotting  coL  =  n  co,,  i.e.  as  multiples  of  co,,  we  would  then  have 
Table  A-l. 
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TABLE  A1 

n 

fin  (0,) 

1 

4.3  x  10n 

2 

6.8  x  1010 

3 

2.7  x  1010 

10 

2.6  x  109 

20 

7.7  x  108 

30 

4.2  x  108 

40 

3.1 x  108 

50 

2.7  x  108 

60 

2.7  x  108 

70 

3.1  x  108 

80 

4.5  x  108 

90 

1.1  x  109 

100(=  N) 

4.6  x  1010 

Prom  this  table  we  see  (a)  there  is  an  order  of  magnitude  drop  from  n  =  1  to  n  =  2  in  fln  ©,). 
(b)  the  integrand  bounces  back  but  is  always  an  order  of  magnitude  (or  two  or  three)  below  the 
initial  value.  We  are  therefore  led  to  compare,  numerically,  the  value  of  INT.  given  above, 
namely  3.7  x  107,  to  the  approximation 
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V ("■) ‘  —*■  r~  „  ,5/3  -I 

©4  +  N  -  coj 


(0, 


(0(  N 
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=  4.3  x  10. 


(A-29) 


This  is  acceptable  agreement  and  thus  we  consider  it  appropriate  to  accept  the  approximation. 
When  ©l  =  *Sc  =  N/q  from 


N  -  to  =  qk  ~  N 

I  X 


(A-30) 


calling  Akx  =  kx  =  N/q  in  (A-27)  and  using  INT  =  y  (©t)  above,  we  obtain 
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(A-31) 


where  the  (A- 29)  bracket  term  was  used.  Thus 


1 

2 


ue  (N)  =  2q2 


(A-32) 


and,  from  (A-15)  ^  u  2  (n)  =  a'  e  N  1  =  (w,  q2)/N  (from  e  =  cot  q2  and  a'  ~ 


1),  the  ratio 


iu2(N) 


(A-33) 


which  is  tlie  same  as  the  frozen  case  (A- 17)  within  constants  of  order  unity. 

This  proves  that,  as  far  as  our  purposes  are  concerned,  the  "frozen  waves"  hypothesis  is 
valid.  Finally  it  will  be  noted  that  our  conclusions  were  never  based  on  the  absolute  value  of 
(A- 17)  but  rather  the  ratio  of  this  quantity  taken  from  the  atmosphere  and  ocean.  For  this 
reason,  plus  the  fact  that  the  constant  of  order  unity  involved  would  be  universal,  the 
conclusion  is  independent  of  the  actual  value  of  that  constant. 
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Notes  Added  in  Proof 


1  -  The  unusual  history  behind  the  theory  of  Kolmogorov  and  Obukhov  regarding  inertial 
range  turbulence  is  given  in  "Statistical  Fluid  Mechanics"  by  A.S.  Monin  and  A.M.  Yaglom, 
MIT  Press,  1965,  vol.  1,  pp.  12-17.  In  addition  to  the  rhyme  of  L.F.  Richardson  mentioned  in 
the  text  Kolmogorov  was  also  aware  of  the  "4/3”  length  scaling  law  of  the  diffusion  by 
turbulence  of  clouds  which,  as  was  later  realized,  follows  dimensionally  from  the  cascade 
hypothesis. 


3  -  The  present  report  does  not  discuss  the  problems  arising  from  variable  N  and  variable 
shears  with  altitude,  their  filtering  effects,  etc.  In  regard  to  saturation  see  Van  Zandt,  T., 
Fritts,  D.,  Pageoph.,  130,  p.  399  (1989)  who  examine  some  of  such  effects.  The  growth  of  wave 
amplitude  with  altitude  is  a  crucial  factor. 


4  -  The  agreement  between  the  model  and  the  e  measurements  of  Hocking  and  yu  (to) 
measurements  of  Balsley  and  Carter  that  is  mentioned  on  p.  21  was  not  used  during  my 
construction  of  the  present  model.  It  was  learnt  by  the  author  only  at  the  last  stage  of  writing. 
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In  other  words  I  claim  that  it  was  really  a  "prediction”  that  was  consistent  with  an  unexpected 
observation.  This  is  the  reason  that  I  think  that  simultaneous  measurements  of  yu  (co)  and  e 

are  both  justified  and  worth  the  effort. 

5  -  Regarding  the  simultaneous  presence  of  cascade  and  saturation  it  should  be  clearly  stated 
that  this  theory  assumes  (a)  cascade  processes  control  co  and  kx  dependent  PSD's,  and 
(b)  saturation  controls  the  kz  PSD  of  horizontal  velocity  fluctuations.  A  planned  future  report 
will  describe  the  simultaneous  presence  of  these  two  processes  in  detail  and  will  explain  why 
(b)  does  not  significantly  interfere  with  (a). 
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